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Abstract 

We study the asymptotic regime of strong electric fields that leads from the Vlasov-Poisson 
system to the Incompressible Euler equations. We also deal with the Vlasov-Poisson-Fokker- 
Planck system which induces dissipative effects. The originality consists in considering a situation 
with a finite total charge confined by a strong external held. In turn, the limiting equation is set 
in a bounded domain, the shape of which is determined by the external conhning potential. The 
analysis extends to the situation where the limiting density is non~homogeneous and where the 
Euler equation is replaced by the Lake Equation, also called Anelastic Equation. 
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1 Introduction 


1.1 The Vlasov-Poisson equation in a confining potential 

We are interested in the behavior as e tends to 0 of the solutions of the following Vlasov equation 


dtfs + • Vai/e — -Vx^ext + — 0, 

£ 


where the potential is defined self-consistently by the Poisson equation 

A^^s = -^Pe, Peit,x) = J fe{t,X,v)dv, 


(V) 

(P) 
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Figure 1: Snapshot of a 2D simulation of confined charged particles. Particles are subjected to the combination 
of a harmonic and isotropic external potential, a strong Coulomb repulsion, a friction and a noise. An external 
force has been added from the left to the right in the lower half of the cloud, in order to set the particles in 
motion. Left: instantaneous locally averaged density field. The density is almost uniform inside a ball, and 
almost zero outside. Right: instantaneous locally averaged velocity field. (By courtesy of A. Olivetti [34]. 1 


and where e ^4>ext is a strong external potential applied to the system. The problem holds in the 
entire space: x G v G and it is completed by an initial data with finite charge 


fe 


t=0 


^init 

J £ ’ 



du dx = m S (0, oo). 


( 1 ) 


Notice that is of size and we shall consider the applied potential j4>ext also of size 
The problem is motivated by the study of non neutral plasmas (see [12] for a review): these are 
collections of particles all with the same sign of charge, for instance pure electron, or pure ion 
plasmas. There are several methods to confine such a plasma, among which the Paul trap, which 
uses an oscillating electric field. The Penning trap, which uses a combination of static electric and 
magnetic fields, is also standard, but (|v| is not directly relevant to this situation since there is no 
magnetic field in it. A non neutral plasma picture has also been used to describe trapped neutral 
atoms [31] . in the regime where multiple diffusion of quasi resonant photons induces an effective 
interaction force between atoms which is formally similar to a Coulomb force [10] ■ In this case, 
the system is however dissipative; a standard way to take this effect into account is to add to 0 
a Fokker-Planck operator acting on velocities [1]. We will also discuss this situation. In these 
physical examples, the small e limit is indeed relevant in many experimental situations. Figure 
corresponds to a numerical simulation of such an experiment. It strongly suggests the existence of a 
limiting fluid model where the density is nothing but the characteristic function of a ball. Our goal 
is to justify that, indeed, a simpler model, purely of hydrodynamic type, can be used to describe 
the particles in this asymptotic limit. 

In fact, we shall see that the limiting model holds in a domain the shape of which depends on 
the external potential d>ext- But, to start with, we can consider a quadratic and isotropic potential, 
say: 

^eM=^\x\^ ( 2 ) 

where we remind the reader that N stands for the space dimension. It corresponds to the case 
displayed in Figure [l] The confining potential e“^$ext tends to strongly localize in space the 
particle density. On the support of the limiting density p, the electric force £“^Va;$ext + 
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should be of order one. By 0. this imposes that A$ext + ■ {s ^Va;<i)ext + Vx'he) = 

Ad)ext — Pe = 1 — Pe is of Order ^(e) on the support of p for the potential (H). Clearly, due to 
the condition of finite charge 0’ the limiting density cannot be constant uniformly on the whole 
space. The intuition is that the limiting density has the same radial symmetries as both the external 
potential ([^ and the Poisson kernel, see (14) below. Actually, we shall prove some convergence of 
Pe to 

ne(a;) = lB(o,_R)(a;), (3) 

where lu denotes the characteristic function of the set U. The radius R depends on the total 
mass m so that the charge constraint Q is fulfilled. In order to find a hydrodynamic description 
of the particles, it is convenient to associate to the particle distribution function the following 
macroscopic quantities 

Current: Je{t,x)'^= J v fe{t, x,v) dv, 

Kinetic pressure: Pe(t, x) J u 0 r:/^(t, x, u) du. 


It turns out that the current looks like 


Je{t,x) = Pe{t,x)Veit,x) -^ ne{x)V{t,x) = lB( 0 M{x)V{t,x), 

£->■0 ' ' 


where V solves the Incompressible Euler system in B{0,R): 


dtV + V, • (K 0 K) + = 0, 

V, • K = 0, 


(4) 


(IE) 


with an appropriate initial condition, and no flux boundary condition on dB(0^R). In (IE), the 


pressure p appears as the Lagrange multiplier associated with the constraint that V is divergence 
free. This incompressibility condition comes from charge conservation: integrating ( [V| with respect 
to the velocity variable u, we get 

dtPe + Vx • de = 0. (5) 

Letting e go to 0, with 0 and 0 , we deduce that V is solenoidal. Obtaining the evolution equation 
for V is more intricate. 

The analysis of such asymptotic problems goes back to [S], where a specific modulated energy 
method was introduced. It has been revisited in [29], still by using a modulated energy method, 
but which is able to account for oscillations present within the system. Accordingly, more general 
initial data can be dealt with in [^. However, these results hold either on the torus T^, or in 
the whole space with data having infinite charge, that is JJ f{x^v)dvdx = oo. A case with finite 
charge, but a different Poisson equation which leads to a compressible hydrodynamic limit, has 
been considered in |19j . again with a modulated energy. Our goal in this article is twofold: 


To prove the convergence to (IE) in the case of a trapped system, with finite charge. Even 
though our proof also relies on a modulated energy functional, there are new difficulties: the 
shape of the domain on which the limiting equation (IE) holds is determined by the external 
potential 4)ext, and a careful treatment of the boundary is needed. 


To prove the convergence to the analog of (IE) in the case of a trapped dissipative system. 


Both improvements are relevant for experiments on non neutral plasmas or large magneto-optical 
traps. 

1.2 Statement of the results 


In what follows we shall deal with a smooth solution (t, x) V (t, x) S (possibly defined on a 


small enough time interval [0, T]) of the incompressible Euler equation (IE) set on the ball B{0, i?). 
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completed with no-flux boundary condition 


V(t,x) ■ vix) = 0, 

\x\—R 


( 6 ) 


where i'{x) denotes the outward unit vector at a; G dB{0, R) (namely i'{x) = a;/|a;|). We work with 
solutions V that belongs to L°°{0,T] R))), for a certain s > 0 large enough. 

Theorem 1.1 r[38|-|39]l Let 1/“'* : B{0,R) —> be a divergence free vector field in iJ®, with 

s > 1 -I- N/2, satisfying the no flux condition ■ i/ = 0 on dB(fl, R). There exists T > 0 and a 
unique solution V G L°°{0,T-, H^{B{0, R))) of (IE) with the no flux condition (|^. Moreover, we 
have 

sup (||E(t)|U. + WdtVmH^-^ + l|V,p(t)||ff. + < C(T) 

0<t<T ^ ' 

for some positive constant C{T) depending on T and the initial datum. 

UN = 1, the only divergence free vector held f/™** satisfying is yinit ^ Q then the solution 
given in Theorem |1.1| is V = 0. 


For further purposes, we need to consider an extension T" of the solution V to (IE) with (k3|), 


defined on the whole space and compactly supported. Namely we require 'f G L°°{0,T-, 
to satisfy 




B{0,R) 


= V, 


r 


\x\—R 


= 0. 


)) 

(7) 


= 0, 'Tit, x) ■ v(x) 

1SL"\B{0,2R) 

For the construction of such an extension, we refer to HZl Chapter I: Theorem 2.1 p. 17 & Theorem 
8.1 p. 42]. For an extension which is in addition divergence-free, see Lemma B.l in the appendix. 

In order to state our first result, we need to introduce an auxiliary potential function $e. 
Suppose that indeed holds true. Then, by using © and A$ext = 1 for the potential § , we 
infer, for £ —>■ 0, 


A(<I>ext + e4)e) — A<I>ext — Pe ^ A$ext — 

Moreover, since we want the electric force £“^V 2 ,<i)ext + = £~^(Va;<i)ext + to be of 

order one on the ball B{0, R), this imposes 4)ext + to be close to a constant, say zero, on the 
ball B{0, R). It is therefore natural to look for a solution $e to the Poisson problem 


A$e(a;) = 1 - ne{x) = 1 mn\b{o,r), $e = 0 in B{0,R). 

In this specific case, we can find an explicit radially symmetric solution: 

R^ 


( 8 ) 




27V 7V(7V-2)|x|^-2 2{N - 2) 

\x\^-R^ R\ „ 

--^fo(k|/7?) 

\i\x\ - Rf 


UN >2, 
UN = 2, 
if TV = 1. 


(9) 


With $e and rze in hand, we split the Poisson equation ^ as follows, where rie is defined in ([^, 

A a: ^ X 1 ~ Reix) nJx) — pJt.x) 1a* 1**/X 1** / X 1a* 

A,$e(C x) = -^ + ■ - - - - A$ext = - A,$e(a:) + x) - - A$ext, 

£ see y £ £ 

namely, we have 

$E(a;)-I--^'ext =-4‘e(x)-I-^4'e(t,x), A„;'^/s{t,x) =-^{ne{x) - pe{t,x)), (10) 

e s y/s y/e 
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where represents the fluctuations of the potential. According to [5], we introduce a modulated 
energy: 

^ fedvdx+^ J |Va;^'epdx+ ^ JJ $e fedvdx. 

When the external potential is given by ([^, we shall establish the following statemenl]^ 


Theorem 1.2 Let R)) satisfy = 0 and the no flux condition ja. Denote 


by V the solution, on [0,r], to (IE) with the no flux condition (|6| 
Y a smooth extension of V satisfying the conditions 0. Letfl 
of integrable functions that satisfy the following reguirements 


given in Theorem 1.1 

,N 


Consider 


X —)■ [0, oo) be a sequence 


{ jj di; dx = m, 

1 // 1^- ^“‘1' ff^dvdx+\ J fr^dvdx^ = 0. 


Then, the associated solution fg of the Vlasov-Poisson equation 0-0 satisfies, as e —>■ 0, 
i) pe converges to in C°([0,T];— weak — *); 
a) converges to 0 uniformly on [0,r]; 

Hi) Jg converges to J in .^^([0,r] x M.^) weakly-*, the limit J lies in L°°(0,T;L^(IR^)) and 
satisfies J\[o,t]xB(o,r) = V, V, ■ J = 0 and J ■ = 0. 


Remark 1.3 (i) Here, we were not very precise about the type of solutions to the Vlasov-Poisson 
system (|v| -© we are considering. We refer to 1351 \28l for the construction of global regular 
solutions to the system and some extra conditions to ensure the propagation of regularity. There 
are also weaker notions of solutions (weak solutions or renormalized solutions) to which our theorem 
can apply. We refer the reader to the introduction of \2SI) for a discussion about these solutions. 

(a) The second part of the hypothesis 0 imposes that the initial modulated energy is small; 
this is a strong hypothesis on the initial data. When the problem is set on the torus, or on the whole 
space with infinite charge, it can he relaxed, see \29f . In the present framework, going beyond 0 
would certainly require a fine description of boundary layers on {|x| = R}. Assuming (11), point 
ii) of the theorem then ensures that the modulated energy remains small at later times. As typical 
initial data satisfying ([II|), we can take 




nfix) — SgAx(x) ^ f V — 

Gr 


tN 


where x € R)) and where G is a nonnegative function that belongs to the Schwartz space 

and satisfies J Gdv = 1 (for instance, G is a normalized Gaussian G{v) = (2Tr)~^/^ exp(— |np/2 )). 
Then, we choose tTg —>■ 0 as e —)■ 0 and Sg = o(y^) (so that — SgAx > 0 for e small enough). In¬ 
deed, we easily obtain ff <i)e dv dx = 0, ff dv dx = ag(f Ivl^G(v) dv)(f riedx)—> 

0 and 'kg = 5g/y/ex, hence f da; = dlje f da; —>■ 0. 


We wish to extend this analysis by dealing with more general external potentials. We distinguish 
two situations depending on the expression of the external potential: 

• The quadratic potential 

N 2 

^ext(a;) = (12) 


^Throughout the paper, we denote by the space of bounded measures on X C R^. It identifies with the dual 

space of the separable space C'o(X) of the continuous functions that vanish at infinity. 
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Figure 2: Some level sets of $ext and the support of rig- 


with Xj > 0, 1 < j < N, in dimension > 2 is typical to model non neutral plasmas [T^ 
or magneto-optical traps experiments. In this case A<i>ext is still a constant, that therefore 
determines the value of the (uniform) particle density Ue on its support. But the problem has 
lost its symmetries and the shape of the support becomes non trivial. We shall see that pe 
tends to a uniform distribution rig, supported in an ellipsoid. However, we point out that the 
support of rig does not coincide with a level set of $ext- An example with N = 2 is given in 
Figure The potential <I>e can be computed rather explicitly, and Theorem |1.2| generalizes 
directly. See Section [2T] for a precise statement. 


• In the case of a non quadratic potential, under suitable hypotheses on <I>ext, the limiting 
density rig still has a compact support /C and is still given on /C by rig = A<I>gxt. However, 
rig is clearly no longer constant on K,. The identification of K, and rig relies on variational 
techniques, with connection to the obstacle problem. It is still possible to prove the analog of 
Theorem 1.2 but, since rig becomes non homogeneous, instead of (IE) the limiting equations 
are now the so-called Lake Equations, see e. g. [50]: 


r atE + t/-V,E + V,p = 0, 

\ V, • (rigE) = 0. 


(LE) 


Such model — also referred to as the Anelastic Equations — arise in the modelling of at¬ 
mospheric flows [35]; we refer the reader to m for the justification of a derivation from the 
compressible Navier-Stokes system. As a matter of fact, we can observe that the first equation 


in (LE) maybe writen in the following conservative form 3t(neE)-fV2,-(ngE®E)-|-rig = 0. 


The construction of $g and JC, and a precise statement of the corresponding convergence the¬ 
orem can be found in Section [2?2l 


Motivated by actual experiments, we will also generalize the results to the case where a Fokker- 
Planck operator acting on velocities is added to Eq. 0 - Our starting point then becomes: 


dtfe + V ■ V,/g - • V„/g = Lf, 


(VFP) 
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with 


Lf = V,- {vf + 0V„/) = 0V„ • Mo,eV, 


/ 


Mo^eiv) = 


o-\^v/m 


MofiJJ" (27r6»)^/2 

for some 0 > 0. Equation ( |VFP ) is still coupled to the Poisson equati on ([P ). Using a modified 
modulated energy, we are able to show in this case that solutions of (VFP) and © with well- 
prepared data converge when e and 9 tends to 0 in the sense of Theorem 1.2 to UeV, where V is 
now the solution of the Lake Equation with friction 


dtV + V- V,U + + U = 0, 

■ {n^V) = 0. 


(13) 


On the boundary, we still have the no-flux condition®. For the sake of completeness, the necessary 
analog of Theorem |1.1| for the systems (LE) and ( |13[) is s ketched in appendix See Section for 
a precise statement on the asymptotic behavior of (VFP) and its proof. 


2 The limit density Uq and total potential $e 

As said above, we have a clear intuition and explicit formulae for the equilibrium distribution rie 
and the potential <i>e in the specific case of the isotropic external potential ([^. Let us discuss in 
further details how <i>ext determines rie and its support, and how the auxiliary potential <i>e, which 


plays a crucial role in the analysis through the decomposition (10), can be dehned. 


We remind the reader the definition of the fundamental solution, hereafter denoted P, of (—A) 
(mind the sign) in the whole space 

1 


-n / \ def 

r(a;) = 


iV(iV-2)|BR«(0,l)|-|a.r-2 
In |x| 


if V > 2, 
if V = 2, 
if V = 1. 


(14) 


2.1 The case of a general quadratic potential 

Let us consider in this section the case of a quadratic potential ( [T^ . We have A$ext = ^ 

which is constant in space. It gives the value of the equilibrium density on its support since we still 
expect pe —>■ l;cA$ext- But it remains to determine this support Supp(ne) = 1C C on which we 
have the volume constraint 

^ 1 

Bdx= |/C|^-2 

coming from ([^. Note that a quick computation reveals that K, can be neither radially symmetric, 
nor a level set of $ext- 


m = 


1.2 


In order to extend Theorem 
and a function —)■ M such that 

N 


for a potential as in (12), we need to construct a domain K, C 


pN 


A$e(a:) = W 1: 


.=1 




<i>p = 0 in K.. 


(15) 


The starting point is the observation that given a = (oi, ...,ajv) G (R+)^, then the characteristic 
function of the ellipsoid 

N 

ICa = {x G 
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generates an electric potential which is quadratic inside the ellipsoid. This can be found for instance 
in Chapter VII, § 6]; the computation there is for iV = 3, but the extension to the case V > 3 is 
straightforward, and the two-dimensional case is treated by using arguments from complex analysis 

in [T7] . 

For X S , we denote by cra(a;) the largest solution of the equation 


N 


E 


a]+^ 


= I 


(with G M as unknown). Consequently, x G K-a holds if and only if aa{x) < 0. By convention, 
1 ^ 0 ( 0 ) = — 00 . This quantity can be seen as an equivalent of the radial coordinate in the ellipsoidal 
coordinate system. It allows us to construct a solution to (15) where JC is an ellipsoid, the coefficients 
of which depend on the mass m and the A^ ’s. 

Proposition 2.1 Let a = (ai,...,aAr) G 


0) 


If N > 3, then 


r — 


N 


X 



- 1/2 


IlK' + s) 


IlK' + s) 

1=1 


ds 


- 1/2 


ds 


if CTaix) > 0, 

if CTaix) < 0. 


(ii) If N = 2, then 


r-klK:^{x) = - (0102)X < 


- In aa{x) -I- ^ ^ + V (a? + tTa(x))(a‘^ + cra{x)) 


f + OO ^ 
J(Tn.{x) A — ^ 




ds 


a(^) ^ ® + s){al + s) 


-In ( + 02)2 j _ 


^+00 2 




ds 


+ s a /{af + s)(a 2 -I- s) 


if aa{x) > 0, 
if o-a(x) < 0. 


Remark 2.2 A n alternative point of view for the two dimensional case is to work with the electric 
field instead of the potential. We refer to EE for expressions of the electric field generated by 
ellipses in N = 2. In the case of a uniform charge distribution, the electric field is linear inside the 
ellipse, with the same coefficients for the quadratic terms as those coming from the expression in 
(ii). 


We define the mapping Z : ^ hy 


r+oo 

2M) = / - 

Jo 



- 1/2 


ds > 0. 


(16) 


From Proposition 2.1 we know that the potential generated by Ijc^ is quadratic inside Ka, up to an 
additive constant. The coefficients of the quadratic terms are the “(nj=i (^j)2k{a)/A, 1< k < N. 
The idea to make the connexion with the external potential $ext is now to adapt the a^’s so 
that the quadratic terms in T (inside JCa) cancel out the quadratic terms of ^ext, so that 

(A$ext)r ★ IjCa. + ^ext is coustaut in JCa- We observe that related to the total charge of 

the ellipsoid Ka since 


N 


N 


N 


m = 


= \iCa\Y. xf = irr»(o, 1)1 n «d E V'- 


i=i 


V/=1 






















We shall thus need to solve equations in a of the form Z(af,a%) = z, where 2 G (IR!j_)^ is given. 
Therefore, we are interested in showing that Z : —>■ is a smooth diffeomorphism. 

When TV = 2, explicit computations may be carried out. 


Proposition 2.3 Assume N = 2. Then, for any a G 


Z{a) = {Z,{a),Z2{a)) 


ai + y/oiof Oil + ^OL\ai 


Moreover, Z : —>■ (]R!j_)^ is a smooth diffeomorphism and its inverse is given by 


Z-\z) = {{Z-^Uz),{Z-^Uz)) 


2zi 


2z\ 


21(21 + 22 )’ 22(21 + 22 ) 


Proof. The explicit formula for Z{a) comes by computing the Abelian integral 


^+00 


ds 


(ai +s)3/2(a2 + s)i/2 




Oi 


02 


ds 


2 

oi + \J0i\a2 ' 


for oi 0 . 2 , and the formula holds true when oi = 02 as well. The formula for the inverse then 
follows by direct substitution. □ 

For N > 3, we no longer have simple expressions for Z. However, we shall prove that Z : 
(K^)^ —)■ (M^)^ is a smooth diffeomorphism by using the fact that Z : —)■ (K^)^ is a 

gradient vector field associated with a strictly concave function. 

Proposition 2.4 Assume N >2 and let us define the function ( : —>■ K by/ 


C(a) = 



ds 


if N >3 
if N = 2. 


Then, C : (M!j_)^ —>■ K is smooth, strictly concave and it satisfies Vf = Z. Furthermore, = 
Z : (K^)'^ —)■ (K^)'^ is a smooth diffeomorphism and for any 2 G (K^)'^, Z~^{z) is the unique 
minimizer for 

inf (2-a-C(a)). (17) 

In ( |17[ ) we recognize the minimization problem that defines the Legendre transform of f. This 
gives a way to compute numerically Z~^(z) through the minimization of a convex function. 


Proof. The smoothness of f is clear and = Z follows from direct computations. If TV = 2, the 
strict concavity of C is straightforward and the fact that V( = Z : (K^)2 ^ (M^)2 is a smooth 
diffeomorphism comes from Proposition 2.3 for any 2 G (ffi+)^. 


for any 2 G (ffi+)^, Z~^(z) is a critical point of the 
strictly convex (since ( is strictly concave) function a i—>- z ■ a — ((a), hence is the unique minimizer 
of that function. We assume now N > 3. Then, for each s G K+, the function 


N 


- 1/2 


G (M!j_)^ I-G ]^(afe 


\k=l 


is logarithmically strictly convex since lnot37s(Q;) = (—1/2) + s) and Hess(ln otUs, a) = 

(l/2)Diag((ai + s)“2^ -|- s)“ 2). Consequently, —((o) = /q^°° z!7s(a) ds is a strictly convex 

function of a. Let us show that the Jacobian determinant of Z never vanishes, that is Hess(/', a) 
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is everywhere negative definite. For that purpose, for v G , we write —ri^Hess(^, a)r; = 
Jg ?;^Hess(t37s, a)?;ds, and thus it suffices to show that Hess(t37s,a) is positive definite for any 
s > 0. Now, we write vusia) = exp(lnot37s(a)), thus d^f.Ws^a) = exp(lnoti7s(a))[9|j,(lno-n7s)(a) + 
9j(lnoci7s)(a)9fc(lnoci7s)(Q;)]. Therefore, if u ^ 0, we obtain 


z;^Hess(-ci7s, a)v = vjs{a) 


N 


u^Hess(lnoti7s, a)u + Vjdj(liiO'cUs){a) 


vi=i 


N 2 


>-!x7s(a)'U^Hess(lnoros,a)'(; = ^ > 0, 


as wished. 

Let us now fix z € and consider the minimization problem In view of the negativity 

of (^, this infimum /i belongs to [0, +oo). Since C is strictly concave, this problem has at most one 
minimizer. Let us show that it has at least one by considering a minimizing sequence (a”)n>o G 
(K^)^. We claim that the sequence (a”)„>o is bounded. Indeed, we have z-a" —C(a") -G n G K+, 
and since C ^ 0; this implies z • a" = C(q^") + I* + o(l) < /r + o(l). Using that all components of z 
are positive, the claim follows. As a consequence, we may assume, up to a subsequence, that there 
exists a = (oi,..., ajv) G such that a” —>■ a as n —>■ +oo. In particular, C(a"') = z ■ + o{l) 

converges. We now prove that at most two components of a vanish. For otherwise, Fatou’s lemma 
would yield 


r+oo 


N 


- 1/2 


+00 = 


]J(afc + s) 


^+oo 


N 


d5 = 


, - 1/2 

liminf ( TT (a^ + s) | ds 

n—>-+oo \ / 


\k^l 

/•+00 / ^ 


\k^l 


- 1/2 


< liminf 

n—>-+oo 


0 




\k=l 


ds = liminf(—^(a")), 

n—>-+oo 


contradicting the convergence of (('(a"))„gN- It remains to show that a has no zero component to 
ensure that [i + o(l) = z • a" — C(q^") z ■ a — C(q!) so that a G is actually a minimizer 

for ( [T7| ). We then assume that ai = 0, for instance, and show that for sufficiently small 5 > 0, 
z- (d, 02 ,..., aAf) — C(<I) ci 2 ,..., Oat) < z • (0, 02 , c^n) — C(0, c« 2 , ciAf)- This reaches a contradiction 
for n large enough. We thus compute 


D{6) = (^z ■ ((5, 02 ,..., oiv) - as, 02 ,..., oiv)) - (^z ■ (0, 02 ,..., oat) - C(0, 02 ,..., oat)) 


/-loo / N \ 

= zi(5 + J (d + s) f (ofe + s) I 


- 1/2 


r+00 / N \ -1/2 

«-i/2 ( TT (rvi^ -\- 1 d.5 


-a-ir. 


0 sl/2(5 + s)l/2[sl/2 + (5 + s)l/2] 

by monotone 
j'+oo (nf=2(“fc + ®)) 

/o sl/2(d + s)l/2[sl/2 + (5 + s)l/2] "" 'Jo 


ds - / s"!/^ I ]^(afc + s) 

\fc=2 2 

- 1/2 

ds'l. 


N 


17 ^ 11^ + ®) 


\k=2 


As (5 —I 0, we have, by monotone convergence, 

,- 1/2 

ds 


-t-00 1 /IV 

2s3/2 


- 1/2 


^(n(afc + s)) ds = +oo. 


k=2 


hence for <5 sufficiently small, D{S) < 0, as claimed, 
for ([T7|). It then follows that VC(o) = z as wished. 


We may now construct a solution to (I5|. 


Therefore, o G and o is a minimizer 

□ 


lO 









Corollary 2.5 (Construction of the function $e for quadratic potentials). Let N > 2 and 

assume that 

«>ext(a;) = 

j^l 3 

with Xj > Oj 1 < j < N. Let us also fix m > 0. Then, there exists a unique a G such that 


N 


\2 - 


m 


and 


m 


1 


2|BR«(0,l)|Er=iA^^ V 1-. 


Zial...,a%)=\^\ .(18) 

k \ 1 / l<i<N 


Therefore, there exists a constant n, depending only on the A,- ’s, N and m such that the function 


N 


— ‘f’ext + I ^ ) f" * + K 


is convex and satisfies 
^ 1 


— A<i)e = I ^ j with, furthermore, $e = 0 m ICa and $e > 0 in \ JCa- (19) 


Proof. We define A > 0 such that A ^ ^ constant k by the formulas 4k = 

_/2 

“A“^ (n^i ®i) (n5=i(®i +'®)) ds if A > 3 and 4k = -A“^(aia 2 ) ln((ai + a2)^/2) if 

N = 2. The existence (and uniqueness) of a satisfying the conditions (18) then ensures that $e = 0 


in ICa and — A<i)e = A Then, from the formulas in Proposition 2.1 we get, in {oa > 0}, 


4A2 


T"r 

11^=1 aj 


■4)e(a;) =r* Ik:,,( a;) + 


4A2 


T"r 

11 ^ = 1 O-J 




4A2 


T"r A^ 

n ^-^1 ^3 


r+oo 

I(Ta (x) 
N 


i-E 


N 2 
^3 


N 


. 1 + 5 

j=i ^ 


- 1/2 

]^(aj+s)l ds 


f* + 00 


N 


- 1/2 


- 1/2 


/ IlK' + s) 


ds 








N 


/O 


(Jaix) ( ^ ^2 _ 

IS'”'""” 


- 1/2 


+00 / ^ 

]^(a^ + s) I ds 
a=i 

ds. (20) 


The last integral is positive if aa(x) > 0 since, when 0 < s < aa(x), ^ 2 ^ — 1 > (x) ~ 

1 = 0. In order to see that is convex, we notice that = 0 in {aa E 0} = JCa and that, from 




(20), we have, when o-aix) > 0, and for any direction a; € S 


N-l 


N 


q2^ ,. I y 




- 1/2 


P(a| + s) 


ds 


ni 


N 

j=i “a 


N 

E 


^3^3 




N 

E 


N 


- 1/2 


W_{a] + aa{.x)) 


U'=i 
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since X;.”=i 


jj=i -ji \'^j + o’a(a^)) = 1, which is indeed >0. □ 

Clearly, the ellipsoid ICa is not a level set of the external potential <i)ext (except when all the 
Aj’s are all equal). It is interesting to study the limiting case of a very asymmetric external 


potential. For instance in iV = 2, we consider a trapping potential (12) with a large aspect ratio 


A = A 1 /A 2 1. Direct computations (using Proposition |2.3[ ) lead to 

Ai 


Oi = 


m 


02 = 


1 + 



A 2 

V 

/l + ^ 


Hence 


O'! _ ^1 _ a2 


02 


Thus, the aspect ratio of the particles’ cloud is much larger than the aspect ratio of the external 
potential: this is an effect of the strong repulsion, see Figure for a typical picture. A similar 
phenomenon occurs in higher dimensions. For iV = 3 with cylindrical symmetry, explicit formulae 
corresponding to our Z function are given for instance in mi- It is easy to check that for a strongly 
oblate external potential (“pancake shape”), the aspect ratio of the cloud is again of the order of 
the square of the aspect ratio of the external potential. We can now state the analog of Theorem 1.2 
for a general quadratic $ext- 


2.5 


Theorem 2.6 Let <I>ext be any quadratic potential (12) to which we associate, by virtue of Corol- 
the ellipsoid ICa and the potential <I>e- Let S H^{]Ca) satisfy Vx • = 0 m ICa 


lary 

and the no flux condition § ondK-a- Denote by V the solution on [0,T] to (|IE[) with the no flux 


condition Q given in Theorem 1.1 and consider a smooth extension of V in satisfying 


the following conditions, where R> Q is such that ICa C B(0,R), 

f'(t, x) ■ v{x) 


y 


= v, 


r 


ICa 


R«\B(0,2fl) 


= 0 , 


dKa 


= 0 . 


Let x —>■ [0,oo) be a sequence of integrable functions that satisfy (11). Then, the 

associated solution f^ of the Vlasov-Poisson equation 0”0 satisfies, as e —>■ 0, 

i) Pe converges to Ue = "^7^) ^K-a in — weak—-k); 

ii) converges to 0 uniformly on [0,r]; 

Hi) Jg converges to J in .^^([0,r] x M.^) weakly-*, the limit J lies in L°°(0,T;L^(M^)) and 
satisfies JI jQ = V, Vx ■ J = 0 and J ■ v{x)\gj^ =0. 

The existence of a smooth extension iF of H on satisfying the above mentioned constraints 


follows from P71 Chapter I: Theorem 2.1 p. 17 & Theorem 8.1 p. 42]. See Lemma [B.lj in the 
appendix for a divergence-free extension. 


2.2 The case of a general potential 


We wish now to extend the above results to a general confining potential. When $ext is not 
quadratic, the equilibrium density cannot be expected to be constant on its support. In turn, 
the limiting equation will be more complicated than the Incompressible Euler system. Besides, the 
determination of the domain A = {$e = 0} is a non trivial issue, and its geometry might be quite 
involved [37]. In the following we write Cl — 1C for the interior of K,. 

The pair (rie. Cl) should be thought of through energetic consideration. As it will be detailed 
below, the total energy of the system 0-0 is 



fl dv dx -I- 


1 


dec -f 


1 

£ 


$ext/3e da;- 
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It is natural to investigate solutions whose energy does not diverge when £ tends to 0. Hence we are 
interested in configurations close to the ground state rig defined by the variational problem where 
only the electrostatic part of the energy is involved: namely, we wish to minimize 

£[P] '= J ^ext(a:) dp{x) + i JJ r{x-y) dp{y) dp{x), 

for a fixed m > 0 over the convex subset made of nonnegative Borel measures p of total 

mass m > 0 such that / $ext dp is finite. This problem, which is often referred to as the generalized 
Gauss variational problem, is quite classical and the basis of the theory dates back to [IS]. We 
refer the reader to [SSI Chapter 1] for the case N = 2, and to [71 Theorem 1.2] when iV > 3 for 
the existence of a minimizer under suitable assumptions on $ext- In what follows, we shall assume 
that $ext fulfils the following requirements: 

hi) ‘hext : t K+ is continuous, nonnegative and satisfies d)ext(a;) —t +oo as |a;| —)■ +oo, 

h2) If = 2 or TV = 1, we have lim ($ext + mr)(x) = +oo. 

|x|—>-+oo 

The following statement collects from IZIIMIIST] the results we shall need. 

Theorem 2.7 We assume that the potential cf>ext satisfies the hypotheses hi) and h2). 

{i) The functional £ is strictly convex on 
(ii) The problem 

mi{£[p] ; (21) 

has a unique minimizer Ue which has a compact support of positive capacity. Moreover, there exists 
a constant C» such that 


( r*ne + $ext > G* quasi everywhere, 
y r*ne + <I>ext = G» quasi everywhere on Supp(ne). 


(Hi) Conversely, assume that po G 


and Go are such that 


( r*po + ®ext > Go quasi everywhere, 

\ r*po + ®ext = Go quasi everywhere on Supp(po)- 


Then, po is the minimizer for (211; pQ = Ug 


( 22 ) 


def 


We then define the potential $e = r*ne + d)ext — G*. The constant G* in (22) is called the 
modified Robin constant and quasi everywhere (q. e.) means up a set of zero capacity (which is a 
bit stronger than to be Lebesgue-negligible); see [371 Definition 2.11]. IfiV = 1, (22) holds point wise. 


Proof. The statements for N > 3 can be found in [71 Theorem 1.2]. When N = 2, we refer the 
reader to [Ml Theorem 1.3 for (ii) and Theorem 3.3 for (iii)]. If = 2, the strict convexity (i) is 
not explicited in [36]. Thus we give proofs of (i) for A^ = 2, and (i)-(iii) for A^ = 1. 

The argument for (i) is that if po, pi G and 6 G (0,1), then 


£•[(1 - 0)po + 0pi] - (1 - 0)£[po\ - e£[p^] 

= ^ JJ T{x - y) d[(l - 9)po + 0pi]{y) d[(l - 9)po + 9pi\{x) 

JJ ^{x-y)dpo{y)dpo{x)-^9 JJ T{x - y) dpi{y) dpi{x) 
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Unless po = Pi, the last integral is shown to be positive if TV > 3 in |3 Lemma 3.1]. The case N = 2 
is dealt with in [361 Lemma 1.8], under the restriction that po ~ Pi has compact support. Actually, 
the method used in [3, which consists in writing r(a;) as an integral of Gaussians can be 

extended to the case N = 2 as we check now. The starting point is the equality (see [H equation 

( 12 )]) 

Therefore, denoting p = po — pi and r = \x — y\ and using the dominated convergence theorem 
(on each of the sets {\x — y\ < 1} and {\x — y\ > 1}), we obtain 


^ JJ dp{y) dp{x) df 


= lim 

T —>-+oo 


= lim / 


T—>-+oo 


1/T 


47r 


^-d?lV2-U'U-y) 


dC dp(?/) dp{x) dt 


lim f /^e-‘l«lV2|/i(^)|2dedt 
-^+oo Ji/T J dlT 


T-f+oo 
1 


27r|^| 


IjT 

■|p(e)i^ de 


where, for the second equality, we use J dp = 0, and for the third one, we write as the 

Fourier transform of a two dimensional Gaussian. This clearly shows that ff r(x — y) dp{y) dp(x) 
is positive unless p = 0, ensuring the strict convexity of £ on ^g]]j(m). When = 1, we argue in 
a similar way by observing that 

/* + 00 -I 

-r= dt. 

Jq 'j2'Kt ^ ' 

Indeed, - 1 = /J' a„(e-“'/2*) du = -/^ (u/t)e-“'/2‘ du, thus 

^ + 00 -I /•+00 -| nT pr /• + 00 


r /*+oo 


2 '/ 2 e drdu= / V^du = ry/^, 


0 ^0 


where we have used the change of variable t = ujy/^. Owing to this relation, we can follow the 
same lines as above: 


JJ\^-y\My) dp{x) = ^hrn^ J ^ Jj _ i) dp{y) dp{x) dt 


= lim 

T—>-+oo 


1/T 


2^/Mi 


= lim / — 


T—y+oo 


= lim 

T—>-+oo 


1/T 

c. 


47r 


47r 


JJ e dp{y) dp{x) dt 

-te/2-^ax-v) d^dpiy) dp(x) dt 

|p(^)P d^dt 






MM M 
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It only remains to prove (ii)-(iii) for iV = 1. This is tackled in [37], but with r(a;) = — In |x|. The 
very same arguments apply to the case r(a;) = — |x|/2. □ 


Remark 2.8 To motivate the above computation, one can remark that for N > 1 and under the 
condition f dp = 0, we have, at least formally, 

JJ ^{x-y)dp{y)dp{x) = y iVA'VPda: = J j^\p{f)\^ d^ 

The minimization of the functional £ is connected to an obstacle problem. This connection is 
explained in details in [571 Section 2.5]. 


Proposition 2.9 If Ue is the minimizer of Theorem 2.1 . then h = T * He is the unique solution to 
the obstacle problem 

To find (j) G IIil,^(M.-^) such that 

J \/(j) • V(g — (/)) dx > 0, 

holds for any g G with g — (j) compactly supported and (p > q. e. 

where 'ip{x) C* — $ext(a;). 

We then define the coincidence set 

K {$e = 0} = {T * ne = a - $ext} 

and claim that 1C is compact. Indeed, as jxj -G +oo, we have — T -kne{x) ^ —mr(x) and $ext(a^) + 
mr(x) S> 1 whatever is the dimension N by hl)-h2), thus /C = {T * rie = C* — $ext} is bounded. 


Moreover, by (22), the set Supp(ne) \ 1C has zero capacity. We give some examples in section 2.3 


below where Supp(ne) C 1C, due to the presence of points or regions where A$ext vanishes. These 
points are precisely defined in m Section 3.6] and called ’shallow points’ and it is shown in this 
paper (see Proposition 3.12 there) that it is possible to pass from Supp(ne) to 1C by simply adding 
these ’shallow points’. This fact is illustrated in section [A3| below. 

For a general potential $ext, the variational viewpoint and the theory of the obstacle problem 
provide a definition for the equilibrium distribution Ue, the domain K, (which is not always the 
support of He) and the potential $e- The regularity of <i)ext is not “transferred” to the solution $e 
or r*ne beyond regularity (see [T3[, [6|) since the Laplacian of these functions is discontinuous. 
In addition, the topology of 1C is difficult to analyse in general: 1C may have empty interior or may 
exhibit cusps. Hence, these regularity issues for both 1C and Ug need to be discussed individually. 
Let us then list the properties, which very likely are far from optimal, that we need to deal with 
the asymptotic regime: there exists s > I + A/2 such that 

HI) K, has a non empty interior Cl and dCl is of class C^. 

H2) $ext G C'®+^(M'^), A$ext is bounded away from zero on 1C. 

The regularity assumption HI) on dCl excludes the presence of cusps in K.. Let us point out 
two regularity results derived from the obstacle problem theory. 

Proposition 2.10 Let $ext be a potential satisfying hi) and h2) and consider Ue the minimizer 


of (21) given by Theorem 2.1 and let K. {<I>e = 0}- 

i) \2,‘lf Assume that HI) and H2) are satisfied. Then Cl = 1C and the boundary dCl is 

ii) f \20f Assume that $ext € Then, T * rie S and rig = ln(A$ext) as 


measures. 
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Remark 2.11 A consequence of ii) is that, in (22), we may replace ’quasi everywhere’ by ’every¬ 
where’ since all the functions involved are continuous. Note that H2) then implies that is 
and bounded from below on JC. 


Remark 2.12 Under the low regularity assumption $ext G and when N > 2, it follows 

from Theorem 2] that F * rie G ■ This prevents the singular part of the measure Ue from being 
a Dirac mass or a finite sum of Dirac masses, since the fundamental solution F is unbounded (if 
N >2) near the origin. This however may happen in dimension 1 (see the examples in section \K^ 
below) and in these cases, the relation Ug = lf 2 (A$ext) (as measures) might not be true. 


Proof. For the first point, notice that 12 = /C by HI). In addition, we also have from [23l Theorem 
1] or m Chapter 2, Theorem 1.1], since we assume that A$ext G and does not vanish 

on /C, that the boundary dD, is automatically of class (and even C'«+i+/3 for any (3 S (0,1)). 
If ‘hext is analytic, then dD, is also analytic. 

For the second statement, we first invoke the regularity result of [13] (see also |6], [14]) saying 
that since $ext G C C'^’^(M'^), then F * rig belongs to (^^’^(K^). Consequently, in the 

distributional sense, the compactly supported measure Ug = — A(F*ne) belongs to L°°(K^). Since 
F * Ue = C* — $ext in H (q. e., hence everywhere by continuity of the functions), we infer that 
Ug = A<l>ext in H. If we make the assumption HI), 912 is of class and we then deduce that 
Xle — (A4*ext)l o as a measure. If assumption HI) is not satisfied, then, as noticed in [201 Theorem 
3.10], it follows from [24l Chapter 2, Lemma A.4] that Ug = A$ext holds almost everywhere in 1C, 
which concludes. □ 


With these assumptions HI) and H2), we can establish the following statement, where we point 
out that the limit problem is the Lake Equation (LEI instead of the mere incompressible Euler 
system, since now the equilibrium distribution Ug is inhomogeneous. We obviously need a smooth 
enough solution to the Lake Equation (LE): we may refer to the works [IH] , [33] (when the domain 
12 is simply connected) and [25] (without simple connectedness assumption on the domain 12), which 
rely on a vorticity formulation a la Yudovitch and are then restricted to the dimension N = 2. We 
provide in the appendix (see Theorem A.l) a well-posedness result analogous to Theorem [1.1] valid 
in any dimension and without simple connectedness assumption on 12. 


Theorem 2.13 Let <I>ext ^6 o potential satisfying hi) and h2) and consider Ug the minimizer of 


def 


(21 1 given by Theorem 2.1 and let K, = {$e = 0}. Assume in addition that HI) and H2) are 
satisfied. Let g H^{fl) satisfy Vx ■ (rieF'”'*) = 0 m 12 and the no flux condition Denote 
by V the solution on [0,T] to the Lak e Equation (LE), with the no flux condition (|^ and initial 
condition given in Theorem A.l and consider a smooth extension ofV™^*' satisfying the 


following conditions, where R > 0 is such that 12 C B{0,R), 


r 


= v, 


r 


M"\B(0,2R) 


= 0 , 


y(t,x) ■ v{x) 


an 


= 0 . 


Let x —>■ [0,oo) be a sequence of integrable functions that satisfy (11). Then, the 

associated solution of the Vlasov-Poisson equation (0-(0 satisfies, as e —>■ 0, 

i) pe converges to Ug in C°(0,T;./#^(IR^) — weak—*); 

ii) M’-y,e converges to 0 uniformly on [0,r]; 

Hi) Jg converges to J in .^^([0,r] x M.^) weakly-*, the limit J lies in L°°(0,T;L^(K^)) and 
satisfies j \~ ■ J = 0 and J ■ v{x)\g^ = 0. 

The existence of a smooth extension oiV follows from [371 Chapter I: Theorem 2.1 p. 17 & 
Theorem 8.1 p. 42]. 

For convex potentials $ext, the only situation where we have been able to check the hypotheses 
HI) and H2) (except the quadratic potentials for which A$ext is constant) is the case of the space 


dimension TV = 1 (see Proposition 2.14) and the case of a radial potential (see Proposition 2.15 
below). 
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2.3 About hypothesis HI) for convex potentials 

For the problem we have in mind, it is natural to assume that the confining potential $ext is smooth 
and convex. In this case, one may think that the coincidence set K, or Supp(ue) is convex. We 
have not been able to find such a result in the literature for a general convex, coercive and smooth 
enough confining potential $ext- Actually, the obstacle problem is, in most cases, set on a bounded 
convex domain G with suitable boundary conditions instead of the whole space . 

For the obstacle problem in bounded convex domains G, we can find a convexity result for the 
coincidence set K, in [m Theorem 6.1] in the specific assumptions that A$ext is constant and with 
the boundary condition T -kUe = 1 + ip = l + C* — $ext on dG. Just after m Theorem 6.1], an 
example is given (in a bounded convex domain G) showing that the assumption $ext smooth and 
strictly convex (and F * Ue > ip on dG) is not sufficient to guarantee that JC is convex. Roughly 
speaking, A$ext is constant for quadratic potentials. 

Turning back to the obstacle problem in the whole space K^, the only convexity result we are 
aware of is [S] Corollary 7], which corresponds to the case where A<I>ext is constant. Extending this 
result to space depending functions A$ext is a delicate issue (see however m Chapter 2, section 
3], which is not sufficient for our situation). 

In the one dimensional case and for a convex potential $ext, there is a simple characterization 
of 1C, as explicited in the following Proposition. 


Proposition 2.14 (The one dimensional case with a convex potential) Assume that N = 
1 and that <I>ext : K —>■ ffi of class , piecewise , nonnegative, convex (i.e. <I>gxt *5 nondecreas¬ 
ing) and that <I>ext(a^) — iti|x|/2 ^ 1 for |a;| ^ 1 (so that hi) and h2) are satisfied). We denote by 
^^6 piecewise continuous function associated with the second order derivative o/<I>ext- Then, 


the minimizer Ue for (211 is given by 


rig = (d$;^t) 


(23) 


where a+ and a_ are defined by the equations 


I = <!>;,,(«+) and -|=$;^^(a_). (24) 

Furthermore, Supp(ne) = Supp(94)g^^) n [a_,a+] and {$g = 0} = [a_,a+]. In addition, the 
potential $e is convex. 


Proof. As a Hrst observation, notice that ( 241 has at least one (possibly non unique) solution since 
<l>g^j is continuous, nondecreasing and tends to > m/2 (resp. < —m/2) in view of our hypothesis. 
If the limit at +oo is m/2, it follows from the convexity of $ext that $ext(a;) — $ext(2/) < m/2, 
cici 

(resp. —oo) (resp. —m/2) Since $g,jj is nondecreasing, and if b+ > a+ also solves m/2 = 
$g,jj(6+), this implies that, on [a+, 6-|_], ^g^j = m/2, thus = 0 and this does not change rig. 


Let us use the characterization (iii) in Theorem 2.7 and look for the measure under the form 


= (94>(,xt) 

if and only if 


]a_,a+[ 


which is piecewise continuous. This function satisfies the mass constraint 


fa+ 


m = 


5®ext dx = 4>/xt(a+) - d’ext(a-)- 


(25) 


Now, let us compute r*ne + 4)gxt in [a_, a+] and investigate under which condition this function 
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Figure 3: The potential $ext and the corresponding measure Ug for example 2 


is constant (in [a_,a+]). Elementary computations give, for a_ < x < a+: 

T-kneix)=-^J dy 

= - 2 d’ext(a+)(a+ -x) + - a_) + - J sgn(y - x)^'g,^^{y) dy 

= - ^d-e^t(a+)(a+ -x)+ - a-) + ^d’ext(a+) + ^$ext(a-) - $ext(a;). 


As a consequence, T * rie + d)ext is constant in [o_, 
Combining this with the mass constraint <!>(,,(.(. (a+) — 
follows that, on [a_,a+], 


a+] if and only if $ext(a+) + d’ext(a-) = 0. 

= m yields the relation (24K It then 


r ★ TIq ~h ^ext — 


def 1 

= 2 (d>ext(a+) + $ext(a-) - a+$e,,t(«+) - a-d'Lt(a-)) 
= i ($ext(a+) + $ext(a-)) - ^(a+ - a_). 


It remains to check that r*ne + $ext > C* in ffi. To see this, note that $e = r*ne + $ext ~ C* is 
convex since its (distributional) second order derivative is equal to the piecewise continuous function 
and $e = 0 on [a_,a+], hence is > 0 everywhere. This finishes the proof. □ 

Let us give some examples illustrating Proposition |2.14| 


Example 1 (ID): If <i)ext is of class and is positive on K, then ng{x) = $"xt(a^)l[a_.o+](a^) 
and is absolutely continuous with respect to the Lebesgue measure. We then have Supp(ne) = 
[a_,a+]. 
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Example 2 (ID): the potential <i>ext is C^, piecewise C^, but is affine on the interval [a-,a+] 
(hence it is not strictly convex), where its slope belongs to ] — m/2, +m/2[ (see figure]^. In addition, 
the second order derivative is discontinuous at a_ and continuous at a+ and is positive 

except on [a_, a+]. In this case, we may still define a± as the unique solutions to $g,(.^(a±) = ±m/2, 
and we have Supp(ne) = [o_, a_] U [a+, o+] C [o_, a+] = {$e = 0} and this is then a disconnected 
set. If the slope in the region [a_, a+] where $ext is affine does not belong to ] — m/2, +m/2[, then 
the support of Ug is an interval as in Example 1. 


Example 1 fits the hypotheses of Theorem |2.13[ but not Example 2 since rig is not bounded 
away from zero (near a+). In particular, for Example 2, we have to face new difficulties in solving 
the Cauchy problem (see Theorem A.ll for the Lake Equation (LE). If in the one dimensional 
situation one can easily check that the support of Ug (instead of K.) is smooth, in a similar higher 
dimensional case, the regularity of Supp (rig) is certainly not easy to analyse since we can not 
rely on the results in [231 Theorem 1] or m Chapter 2, Theorem 1.1]. All these issues motivate 
hypothesis H2). 

Let us give now examples which do not fit the regularity hypotheses required in Proposition 


2.14 These expressions are justified through the characterization (hi) in Theorem 2.7 and simple 


computation of T * rig 


Example 3 (ID): Take the potential $gxt(a:) = |a:|. Then, hypothesis h2) exactly means m < 1. 
In that case, we have rig = m(5o and Supp(ng) = {0} = {dig = 0}. 


Example 4 (ID): Take two reals a < b and a convex potential $gxt which is affine on ] — oo, a], 
on [a, b] and on [6, +oo[. Assume also that h2) is satisfied, that is m < min(d)gx^(+oo), —d>gxt(—oo)). 


Then, rig = min 



t(6) digxt(^i) 

b — a 


, m 5a+min 


m — 


dlext(^) d?gxt(a) 
b — a 


,m 


As a consequence: 

- if —m < < m, then Supp(ng) = {a, b} and {d)g = 0} = [a, 6]; 

- if — < —m, then rig = mi5{, and Supp(ng) = {6} = {d>g = 0}; 

b — a 

■ r ^ext(^) dlext(ll) e jc/Afl nl 

- if - - - > m, then rig = mOa and Supp(ne) = {a| = {d)g = 0|. 


Example 5 (ID): Consider the potential d)ext(a;) = |a:| + x'^/2 + max(a: — 1, 0): 

- if m < 2, then rig = m<5o, P *ng(a;) = — m|a::|/2 and Supp(ng) = {0} = {$g = 0}; 

- if 2 < m < 4, then rig = 2do + l[-m/ 2 +i,m/ 2 -i] and Supp(ng) = [-m/2 + 1, m/2 - 1] = {d>e = 0}; 

- if 4 < m < 6, then rig = 26o + (m/2 — 2)(5i + l[_m/ 2 +i,i] and Supp(ne) = [—m/2 + 1,1] = {d)g = 0}; 

- if m > 6, then rig = 25o + + l[_m/ 2 +i,m/ 2 - 2 ] and Supp(ng) = [-m/2 + 1, m/2 - 2] = {d>g = 0}. 


Examples 3, 4 and 5 show that the single convexity hypothesis on $gxt does not guarantee that 
rig is a restriction of the nonnegative measure 9^d>ext (in the distributional sense). It appears in 
these examples that rig is nondecreasing with respect to the mass m, and thus that we always have 
ng < 9^$gxt in the distributional sense. It is an open problem to determine whether this holds true 
in higher dimensions. Here again, these issues motivate the regularity assumptions on $gxt in H2). 

The other situation where we may verify hypothesis HI) is the radial case (see [7l Corollary 
1.4] for a related result in dimension N > 3 for potentials diext)- Let ipg^^t : K_|_ K_|_ be a 

nondecreasing function of class and piecewise C^. Consider now the potential $gxt : t R 

given by $gxt(a:) = +ext(|a^|)- It is then clear that (/jgxt is convex if and only if digxt is convex. 


Proposition 2.15 (The radial case with a convex potential) Assume that N > 2 and that 
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( 26 ) 


$ext : —>■ K is as above. Then, the minimizer Ue for @ is given by 

ne{x) = lB(o,fl)(a;)A$ext(a;), 


where R is defined by the equation 

m= / A$ext(a^) dx or, equivalently, A|i?(0, (27) 

Jb{o,r) 

Furthermore, Supp(ne) = B(0,R) \ B{0, Rmin), where i?min = 0} < -R. In addition, 

the potential $e is convex. 

Proof. The existence of R is clear. We may have non uniqueness only in the case where $ext is 
constant on a ball B{0,Rq) (of positive radius), since is nondecreasing. The potential $e may 
be searched for under the form of a radial function, and we find the expressions 

$e(a;) = (</5ext(R) - (/5ext(|a;|) + (/7g^t(7?)^(R))lB(o.i^) + (Pgxt(77)r(a;)lB«\B(o,fl), 

where r(R) stands for r(?/) for any y G dB{0, R). □ 

Let us give some examples illustrating Proposition |2.15| 

Example 1 (radial): If</5gxt is of class and (pg^t is positive on IR+, thenne(a^) = lB(o,fl)( 2 :)^‘i’ext(a;) 
and is absolutely continuous with respect to the Lebesgue measure. We then have Supp(ne) = 
B{0,R). 

Example 2 (radial): The potential (pext is C^, piecewise but is constant on the interval 
[0, i?o] (hence it is not strictly convex). It does not matter whether the second order derivative of 
(fiext is continuous or not at Rq. We define i? > i?o > 0 by the relation m = A$ext dx, 

or, equivalently, A|B(_0,1)|R^-Vext(7?) = m. Then, Ue = lB(o,i?)\B(o.flo)^®ext, Supp(ne) = 
B{0,R) \ B{0,Rq) C B{0,R) = {$e = 0} and this set is then neither starshaped nor simply con¬ 
nected. Here again, if (/?ext € C'^, this potential does not fit hypothesis H2) since A<I>ext is not 
bounded away from 0 near Rq. 


Let us give now examples which do not fit the regularity hypotheses required in Proposition 


2.14 These expressions are justified through the characterization (iii) in Theorem 2.7 and simple 


computation of T * rie 


Example 3 (radial): Take the potential (fiextix) = r, that is $ext(a:) = |a;|. Then, A$ext = 
{N — l)/r > 0, He = (N — l)|a;|“^ls(o r), with N\B{0, 1)|R^“^ = m, and Supp(ne) = i?(0, R) = 
{<i>e = 0}. 

Example 5 (radial): Consider the potential (pext{x) = r + max(r — 1, 0): 

- if m < A|i?(0,1)1, then Ug = (fV— l)|a;|“^lB(o r), with R = (m/fV|R(0, l)|)^/^“^ and Supp(ne) = 

.B(0,R) = {$e = 0}; 

- if A|R(0,1)J < m < 2iV|R(0, 1)|, then Ug = (N - l)|a:|“^lB(o,i) -b (m - iV|R(0, l)|)(5aB(o,i) and 
Supp(ne) = B{0, 1) = {$e = 0}; 

- if m > 2iV|R(0,1)1, thenrie = (A-l)|xr4B(o,i)-bfV|R(0, l)|(5^B(o,l)+2(^^-l)|a;^4B(oJ^)\B(o.l), 
where R > 1 is such that 2iV|R(0,1)|(R^“^ ~ 1) + A^|7?(0, 1)| = m, and Supp(ng) = R(0,R) = 
{<i>e = 0}. 


Since we assume (fg^t convex and with 0 as a minimum point, it follows that ifg^t has a right- 
derivative at 0, hence the singularity in l/|a;| at the origin for Ug is the worst we can have. The 
radial Example 5 also shows that we may have Dirac masses on a sphere (of positive radius). 
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Our next results guarantees that K, has non empty interior when the confining potential $ext is 
and convex. 

Proposition 2.16 We assume that 0 is a minimum point o/<i>ext and that the potential $ext is of 
class and convex. Then, there exists rg > 0 such that Brg{0) C 1C. In particular, K, has non 
empty interior. 


Proof. We follow the argument of [241 Chapter 5, Theorem 6.2], where we work on r*ne and 
shall use that it is a solution to the obstacle problem given in Proposition |2.9| with the obstacle 
d^ext- 

We first consider the case N > 3 and notice that Supp(ne) has a positive capacity: we fix some 
a S Supp(ne) such that C* = h{a) + $ext(a) (see (22)). Now, since TV > 3, we observe that (with 
cn > 0) h{a) = r -kne{a) = cn\ ■ *ne > 0 and that 0 is actually a global minimum point of 
$ext, thus C* > <i)ext(a) > d>ext(0) and it follows that ip{0) = C* — $ext(0) > 0. On the other hand, 
h{x) ^ mr(a;) tends to 0 < '0(0) at infinity, thus there exists an i?o > 0 such that h{x) < tp{0)/2 

when |a:| > Rq. For xq that will be close to 0, we let v{x) ipixo) + {x — xq) ■ 'C'ip{xo) be the 
affine tangent to aX xq. Since 0 is concave ($ext is convex), we have 0 < z; in . Furthermore, 
if xq is sufficiently close to 0 (depending on Rq), then V0(a:o) is small (since 0 is and achieves 
a minimum at 0) and thus v > 0(O)/2 > 0 on 55(0, i?o)- Since Av = 0, we may now apply [Ml 
Chapter 4, Theorem 8.3] to infer h < v in 5(0, Rq) (this is a maximum type principle proved using 

the comparison function g min(/i, n)lB(o,i{o) + ^lR«\s(o,flo) formulation of the obstacle 

problem given in Proposition 2.9). In particular, 0(a:o) < h{xo) < v(xo) = 'iP(xq), which means 
that, as wished, Xq G 1C. 

Let us now turn to the dimensions N = 2 and N = 1. Then, it may happen that 0(0) < 0, but 
since h(x) ~ mr(a;) tends to —oo < 0(0) at infinity, the previous argument still applies, n 


If one is able to prove that 1C is convex and assuming that 4)ext satisfies H2), then HI) is 
automatically true. Indeed, any point of dIC has then a positive density and we may then apply 
the regularity result of L. Caffarelli (see e.g., m Chapter 2, Theorem 3.10]) which ensures that 
dCl is of class (hence (7^+^ by H2)). 

We conclude with a result from uni Theorem 3.24] on the topology of 1C valid only in space 
dimension two (the proof uses complex analysis). 


Proposition 2.17 ( [20j ]) We assume N = 2. Suppose that $ext is of class and that its Hessian 
is everywhere positive definite. Then, supp(ne) is simply connected, and equal to the closure of its 
interior. Moreover, if 4>ext is (7^’“ for some a €]0,1[, then dlC is a Jordan curve, for some 
/3g]0,1[. 

The above result does not prevent cusps in dIC, but just says that the boundary dK. possesses 
a parametrization. 


3 Asymptotic analysis 


This section is devoted to the analysis of the asymptotic regime e —)■ 0. We shall point out the 
difficulties and necessary adaptations between the case of quadratic potentials. Theorem |1.2| and 


Theorem 2.6, and the general case. Theorem 2.13 For the existence theory of the Vlasov-Poisson 
equation, we refer the reader to [1] for weak solutions and more recently to [281135j where strong 
solutions and regularity issues are discussed. Further details and references can be found in the 
survey |T5] . 
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3.1 A useful estimate on 


Before we turn to the analysis of the asymptotic regime £ —^ 0, it is convenient to set up an estimate 
that describes the behavior of <i>e close to the neighborhood of dJC. In the isotropic case, <i>ext being 
given by ([^, the potential <i>e is defined by ([^, and we observe that there exists C > 0 such that 


0< (|x|-i?) |V,$e(a;)| < C$e(a:) 


(28) 


holds for any x with |x| > R. More generally, for a quadratic potential (121, we can establish the 
following property, based on the formulas in Section [ZT| 


Lemma 3.1 Let <I>e be the quadratic potential defined as in Corollary 2.5 Let 'C : —>■ be 

smooth, compactly supported and such that 'C ■ i'\q^ = 0. Then, there exists a positive constant C, 


depending only on N, $ext (^nd 'f' such that we have, for any x G 

ir-v^,(x)i<c^,(x). 


(29) 


Proof. Since Y is compactly supported and $e is positive in {aa > 0}, we just need to prove the 

inequality for x close to dICa, that is for cra(x) small. We still define A > 0 so that 

From (20), and by Taylor expansion of the integral, we infer that for 0 < aa{x) ^ 1 and 1 < fc < TV, 

A29fc$e(x) = y (aa(x)-^ +C>(o'a(a^)) I =^^^^^^+0{al{x)). 


U'=i 




Let (x) stands for the vector with components Xfe/a^. In particular, for 0 < aa{x) ^ 1, we get 

V<I>e(a^) {x' 


|V$e(x)|=0(aa(x)) 


and 


|V<l>e(:^^)| \^{x)\ 


+ e>((Ta(x)), 


(30) 


where the unit vector field x i—>■ is smooth near dICa and is the (outward) normal on dlCa. 

Now, observe that 


N 


o = afc ^ 


■ 1 




= 2 - 


Xk 


N 


7aix) 


- E 


j=l {aj+<Xa{x)y 


dk(Xa{x) 


^ X? 


= 2^ - dk(Ta{x) X! T + 0{aa{x)). 

“fe \,=1 s- ) 

Therefore, for 0 < CTq(x) 1 and 1 < fc < A^, we have 


( N 2 \ / / N 2 

^ ^ 1 + 0{al{x)) = icrl{x) f ^ ^ 



Oialix)). 


As a consequence. 


^ X? 


A2$e(a;) = -(tI{x) XI T > 


■ 1 “7 


rlix) 

c 


(31) 


holds for some C > 0. Going back to (30), we arrive at 


r(x) • V$e(a:) ='y{x) ■ X |V$e(a;)| = fix) ■ + e>(cra(a:))^ x 0{aa{x)) 

= {0{<Jaix)) + O(aaix))) X 0{aa{x)) = 0{al{x)) = 0{<^eix)), 
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by (31) and since ■ |-^ vanishes when CTq = 0 in view of the no flux condition satisfied by 


This finishes the proof. □ 


In the more general setting considered in Theorem |2.13[ the result is the following and simply 
relies on the use of a local chart. 


Lemma 3.2 We assume that dfl is of class and that H2) is satisfied. Then, there exists a 
constant C such that, for any x G 


ir-v^,(x)i<c^e(x). 


(32) 


Proof. We have already seen that dH. is actually of class (7®+^. Since <i>e is positive in \ K. and 


y has compact support, by a compactness argument, it suffices to show that (32) holds near any 


point a G dfl. Possibly translating and rotating the axis, we assume a = 0 and that the inward 
normal to at a = 0 is ei = (1,0, ...,0). We let Xi = 0(x_l), where Xx = {x2, ...,xn), be a 
parametrization of dVl near 0, with V0(O) = 0, hence 0(a;j_) = 0{\x±_\^). 

We now consider the function ip : — >• M defined by ip{y) $ 6 ( 2/1 + 0(2/j_),2/j_), where 

yi. = ( 2 / 2 , •■•,2/iv) £ Then, ip{y) = 0 when yi > 0, hence, for 2 < j and 1 < A: < and 

if 2/1 = 0, dkipijj) = dlkPiy) = 0) moreover, df i(p{0,y±) = A$ext(0, 2/j.) in view of the equality 

A$ext(a::) = A,c$e(ai) = (A^t/j-(A_L0)ai(/J + X)^2(^j®)^^ij'^)(2^i in {xi < 0(a:_L)}. 

It follows from these relations that, by the Taylor formula and by using A$ext(0) > 0 and 
yi=xi- 0(a:_L) < 0, 


T{y) = T{y) - <p( 0 , 2 /u) - 2 / 1^1 </5(0,2/±) =yi - t)dl(p{tyi,y±) dt > 


and we deduce 


$e(a:) > 


{xi - 0(a:L))^ 


C 


(33) 


Still by the Taylor formula, we have, for 2 < j < N, 


dMv) = 2/1 / (1 - t)dti,3T{tyi^y^) df = 0{yl) 

JO 


and 


dipiy) = yidlpif^.yi.) + yl [ {1 - t)df'p{tyi,y±) dt = 2/1 A$ext(0,2/u) + 0(1/1). 

Jo 

Now, we write cli$e(a;) = di(p{y) (with y = {xi — 0(a;_L),xj,)) and V_L$e(a;) = ±>p{y) — 

diip{y)V±e{y±), thus 

y(x)-V$e(a;) = fTi(x)5i$e(ai)+7j_(ai)-V_L$e(ai) = yi{x)diip{y)+y±{x)-V ±ip{y)-diip{y)'f'±{x)-V ±Q{y±). 

Note that '^xipiy) = 0{yf). Furthermore, since y ■ = 0 on dfl = {xi = 0(a;_L)} and iy{x) = 

(1, —V_l0(x_l))/|( 1, —V_L0(ai_L))|, we deduce 

y{x) ■ v$e(x) =o{yf) + dMy){[nix) - y±{x) ■ Wxe{x±)] - [yi{Q{xx),xx) - y±{eix±),xx) ■ ^±q{x±)] 
= 0{yl) + 0(|i/i|) X 0(|xi - e{xx)\) = 0(1/?). 

We conclude by using ([^. □ 
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3.2 Basic a priori estimates 

Now that we have in hand the limiting density profile Ug and the associated potential held $e, we 
derive some basic a priori estimates from Q-(|P). 

Using the splitting of Poisson equation as in (10), 0 recasts as 


dtfe + V ■ Va;/e - -Vaj^'e ' V„/e- e ■ \7yfe = 0. 

£ V£ 

Let us compute the time variation of the following energies: 

• Kinetic energy 

d 
dt 

• Leading order potential energy 

d 
dt 


^ fedvdx = jj fgdvdx-^ j j V ■ e fe^vdx, 


fgdvdx = V ■ Va,$e fe du dx, 


• Fluctuations potential energy 


d 1 
^2 


dx = / dx = - hi’. 






da; 


vfe du^ dx = -^ 11^ fe du da;. 


By summing these relations, we conclude with the following claim (which applies for all three cases 
for $ext)- 

Proposition 3.3 The solution Ue,^e = ^^>6 + of 0-@ satisfies the following energy 

conservation eguality 

^ /edudx + <^e fedvdx J | da;| = 0. 

Furthermore, the total charge is conserved 

fe{t,x,v)dvdx = // /e(0, a;, u) du da; = m. 


3.3 Convergence of the density and the current 

We assume a uniform bound on the energy at the initial time, namely 

of <1 /[ dvdx+-^JJ<^g dvdx+^J ‘p da; < oo, (34) 


where solves the Poisson equation (101. Then, Proposition 3.3 ensures that the energy remains 


uniformly bounded for positive times. Thus, possibly at the price of extracting subsequences, we 
can suppose that 

fe^f weakly-* in .^^([0, T] x x R^), = J fedv ^ p weakly-* in ^^{[0, T] x K^). 

Going back to the Poisson equation, we observe that 

Pe — ‘ (^aj'Le) 


where, by Proposition 3.3 Va;4>£ is bounded in L°°(0, T; L^(IR'^)). Consequently, we establish the 
following claim. 
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Lemma 3.4 The sequence pe converges to Ug = p strongly in L°°{0,T; H~^and weakly-k in 
.^^([0, T] X The limit f is supported in [Q,T^]xVl. The sequence J’^ = J vfg dv is bounded 

in L'^(0, T; it admits a subsequence which converges, say weakly—k in ^^{[0,T] x K^); 

the limit J is divergence free, supported in [0, T] x and may be written f vfdv = n^W for some 
w e ^^{[ 0 ,T] X R^). 

Proof. Proposition |3.3| tells us that \v\‘^fs is bounded in L°°(0, T; x Hence, by using 

Cauchy-Schwarz’ inequality, we get 

J \Je\dx < JJ \v\ y/Te \/Te<ivdx < fgdvdx^ J fg dv dx^^ , (35) 

which leads to the asserted uniform estimate on the current. We can thus also assume Jg ^ J 
weakly—* in ^^([0, T] x Furthermore, since the second order moment in v of is uniformly 
bounded, we check that 



Note that p^ and satisfy if- Letting e go to 0 yields 

dtp ^ X ‘ J — d = dtng ^ X ‘ J — d ^ X ‘ J ~ d. 


Thus, J is divergence-free. Finally, since lim| 3 ,|_,,oo $ 6 ( 2 ;) = +00 and the second order moment in v 
of fe is uniformly bounded, {/e, e > 0} is tight, and we can write 



fe di; dx dt = mT = 

i-T 


£—^0 


/o 


Pe dx dt 


/ du dx dt = mT = / pdxdt = T Ugdx 


f du dx dt - 


>0 J an 

rT r 


f dx dx dt 


0 JJR^\n 


He dx dt = 


pdx dt 


/o JQ 

i-T 


/O VQ 
rT 


f dx dx dt = 


He dx dt 


10 J Jn 

rT r 


lo Jn 


pdxdt = 


to Jn 


f dx dx dt. 


to J Jn 


It proves that supp(/) C [0,T] x H x R^, and thus supp(J) C [0,r] x ft. In particular, we note 
that /([O, T]xdnx R^) = 0, and J([0, T] x dVl) =0. □ 

In order to define the normal trace of J over dil (that is the sphere dB{0, R) in the case 
we shall use the theory introduced in [5]. As a consequence of the discussion above, we start by 
observing that J belongs to the set of extended divergence-measure fields over R^, 

see 0 Definition 1.1]. Therefore, according to [HI Theorem 3.1], J admits a normal trace J • 
defined as a continuous linear functional over Lip{'-f, dfl), 7 > 1 (see [H Equation (2.1)]) with 

Jn Jn 


where the function (j) € Lip(^, D) in the right-hand side is an extension of <() G Lipipf, dVL). However, 
by Vx • J = 0 and the support property on J, we can rewrite 


■ ^lan’ ‘^) = 0 + / J ■ = -{^x ■ JA) = 0. 
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Another way to see this is to observe that the normal trace from must be the same as the normal 
trace from \ ft, which is clearly zero since J has support in ft. Consequently, 

J-i'\^ = 0 in [C(0,r;Lzp(7,n))]A 

Remark that this is not a pointwise relation. In particular, it may happen that “ 

is nonzero, but this does not prevent the time integral of J • to vanish. 


3.4 Passing to the limit: modulated energy 

We now study the modulated energy 

^ /edt;da;+ ^ j^ JJ <^e fedvdx, 

where all the terms integrated are nonnegative. Let us compute as follows 


“ dt 


-Y -v^fsdvdx = -^ I - Y ■ da; 


\Y\^ 


by using Proposition |3.3| We thus have 

= J {PeY - Je) ■ dtYdx + J ^fdtPe dx- j Y- dtJe dx. 

Here, we are assuming that the solution of the Vlasov-Poisson system 0-0 is regular 
enough so that we can perform all the calculations that follow. Integrating the Vlasov equation, 
we obtain 

dtJe + Va: • Pg H-H- x^e = 0, (36) 

y/e £ 

where we rewrite 

Q _ TX Tl Ti 

V£ Ve V£ v£ 


and 


Aa,«'eVa;4'e = Va, • (Va^^-e O Va^^-e) - ' 


Combining these relations to the charge conservation <§ and integration by parts, we arrive at 
= j{PeY - Je) ' BtY dx + J Y-Vx dx 


- J DxY : (Pe - Va^^-E 0 Va.^'e) dx + j '^x-Y 

+ - / PeV • Va:<i>edx + / • Vaj^'c dx. 




dx 


Ve 

where DxY stands for the jacobian matrix of the vector field Y. For the last integral, since Ue is 
supported in H, we write it as 

[ v,^edx = o 

Jn Ve 

by integration by parts and using that Va, • {p-YY) = 0 in H and the no-flux condition (|^. 
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Let us set 


A direct substitution leads to 


F-r,e y (u - r) (g) (U - r) /e du = Pe - r (g) Je - Je 0 r + (g Y. 
an leads to 

= J{pey-Je)-{dtr + {r-w,)r)dx 

-/D,^:(P,,.-VA®VA)<i. 

dx+^ J Pe't ■ Va;$e dx. 


( 37 ) 


+ / V, • r 




We shall use the shorthand notation A< B when the inequality A < CB holds for some constant 
C > 0, the value of which might vary from a line to another. As a matter of fact, we can dominate 
the second and third integrals of the right-hand side by 

P.^lloo l^jj\v-r\^f,dvdx + j\v,d>,\^dx^ < WD^rWo. Afr,e- 

Let us distinguish the case of the isotropic potential in order to point out the difhculties. When 
$ext is given by we remind the reader that $e is supported in {|a;| > i?}, radially symmetric 
and increasing in |a;|, see Combining this with ( |2^ allows us to estimate the last term in (|3^ 
as follows: 


1 f 


- / Pef^ ■ Va;$edx 

= 

e J 



1 


e J{\a:\>R} 


• 3^ /12) I 

Pe , , (|x|-i?)|V,$e|dx 


< i 


£ J{\x\>R) 


Pe^e dx 


x| — i? 

iC • x/|x| 


|x| - R 




where we have used that belongs to L°°(R.^) since Y is smooth, compactly supported, and 

iC • = 0 on dB{0,R). For a quadratic external potential (12), we can proceed similarly by using 


Lemma [3.1 1 When dealing with a general potential, we made hypothesis H2) so that Lemma 3.2 
applies and (29) allows us to estimate 


J p^y ■ Va;$e dx J /9ed>e dx < =^r,e- 


Therefore, we obtain 


where we have set 


dt 


,e + 


(38) 


r. =' y(p.-r - J,) • {dtf + (r • V,)r) dx. 


The Gronwall lemma yields 


-^,e(i) < + y ® ‘^'^Xe(r)dT 


for a certain constant C > 0. The assumption 0 on the initial data is that lime_j.o c^^,e(0) = 0. 
Hence, we are left with the task of proving that f* re(r) dr tends to 0 as e —>■ 0. We have 


;(t) dr 


€—^0 


y Jin^r - J) ■ (dtr + {r ■ v,)r) dx dr 

[ f {n,r - J) ■ {dtr + (r • v,)r) dx dr 

Jo Jfl 

{fif.'y — J) ■ VxP dx dr = 0, 


0 JQ. 
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since riei^ and J are divergence free on and their normal trace vanish. □ 


It is worth pointing out that the regularity assumption of the sequence of solutions was only 


made to justify the computations leading to (381. If one consider less regular solutions, we have 


to assume that these solutions were constructed through a regularization procedure and that the 


previous calculations were done on these regularizations and hence (38) will still hold. 


3.5 Identification of the limit 


Let us observe that if the initial datum satisfies (0, then ( [34| ) holds true. Let us first justify i): 
we shall show that f Pcxdx —>■ f Uexdx uniformly on [0,T] as e —?> 0 for any y G We 

start by observing that 


< m 


Pe(i,x)x(x)dx 

holds for any x G C'q(M'^). Next, consider x G The charge conservation © yields 

J Pe{t,x)x{x) = y dtPeXdx = - J Vj; ■ J^xdx = J Je-VxXdx, 


(39) 


hence the uniform bound (35) on implies a uniform bound on f pg(t, x)x(x) dx for 0 <t <T. 


By virtue of the Ascoli-Arzela theorem, the set {t >->■ f Pe{t, x)x{x) dx, e > O} is therefore relatively 
compact in (^([O, T]) for any fixed x € C,1(R'^). This property extends to any x & Cq (M^) by virtue 
of (39). Indeed, for any d > 0, we can pick xs G such that ||x — Xilloo < 5/m. It follows 


that 


Pe{t,x)x{x) = J Peit,x){x- Xs){x)dx + J Pe{t, x)xs{x) dx 


where, owing to (39), the former integral is uniformly dominated by 5 and the latter lies in a 
compact set of C([0,r]). Therefore {t i-G / Ps{t, x)x{x) dx, e > 0} can be covered by a finite 
number of balls with radius 26 in (^([O, T]). Finally, since Cg (R-^) is separable, we apply a diagonal 
argument to extract a subsequence such that f pe{t,x)x{x) da; converges uniformly in C([0,T]) for 
any element x of a numerable dense set in C'o(R^). By uniqueness of the limit, we find 


lim J pi:{t,x)x{x) da; = J Uexdx. 


Going back to (39), we check that the convergence holds for any x € Cq(R'^). 


The manipulations detailed in the previous Section prove ii). In order to establish iii), it is 
convenient to introduce the following functional: given A a non negative bounded measure on 
[0,T] X R^, and p a vector valued bounded measure on [0,T] x R^, we set 


sup • 0 - ^ y A|0p 


where the supremum is taken over continuous functions 0 : [0,T] x 
Prop. 3.4], we have: 


^ According to [U 


Lemma 3.5 (ID) If p is absolutely continuous with respect to A, denoting byY the Radon-Nikodym 
derivative of p with respect to A, we have 


:de{\p) = \ J aivi^g [0,oo], 


otherwise .DY (A, /a) = +oo. 
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Clearly (A, /i) i—>■ (A, /r) is a convex and lower semi-continuous (for the weak-* convergence) 
functional. Let r] : [0,T] —> [0, oo) be a continuous non negative function. Reasoning as in [^, we 
show that J € L°°{0,T; since 


■^{VPe.VJe) = 


< 


1 

2 

1 

2 

1 

2 


lo Jm." Pe{t:X) 

pT r 1 


r]{t) da; dt 


lo Jr" PeitjX) 

pT 


Ir" 


vy/feit,x,v) ^/fs{t,x,v)dv 


B«xR~ 


\vffe{t,x,v) ri{t) dv dx dt 


< 


r]dt 


2 

ri{t) da; dt 


becomes, as e tends to 0 

J(r{r]ne,r]J) < \\v\\lHo,t)- 

Reasoning the same way, we get 


^iPe,Je - Pe'^) = 


\Je-Pey\^ 


dx dt 


Pe 




[ If \v — y\^fedvdxdt < f ^y^gdt. 

Jo JJr"xR" Jo 


It follows that (rie, J — UglC) = 0, which identifies the limit J and ends the proof of iii). 

Finally, we can check that the initial data for the limit equation is meaningful by establishing 
some time-compactness on the sequence Jg. Let 


Wii = {e:[0,T]x 


j,N 


5 AT 


, 0 of class C^, supp(0) c [0,T] x fl, (ne0) = O}, 


which is a closed subspace of the Banach space (endowed with the sup norm for the function 
and its first order derivatives). Multiplying (36) by a function in Wr, we shall get rid of the stiff 
terms. Indeed, for such a trial function 0, we deduce from 

d 


dt 


j Je-Qdx = j Js- 9(0 dx - j Q ■ (Voo • Pe) da; - ^ J Pe& ■ dx. 


(40) 


since 0 • Va;<i)e = 0 pointwise in view of the supports. By using the estimates deduced from 
Proposition 3.3 we observe that the first two terms are bounded in L°°{0,T). For the last one, we 
use the Poisson equation (10) and integration by parts to infer 


1 


1 


/ Pe0 • Vx'I'e dx=—^ UeQ ■ Ve dx - Ax'dt^x^e dx 


V~e 


V~e 


= 0 -b y • V,e(0 • V^r^-e) dx 


= e-^A\^.'i'e\^)dx+ Y. 

l<j,k<N' 


^Xn ^^k^xi. Ike 


where we have used that neO(t, •) is divergence free. For quadratic external potentials, an inte¬ 
gration by parts shows that the first integral is zero (since rieO is divergence free). In any cases. 


the right hand side can be dominated by |!V0||oo||V'I'e||icx,(o,T;L 2 ( 
L°°{0,T). Reporting this into (40) allows us to conclude that 


and it is thus bounded in 


d 

dt 


Je • 0 dx is bounded in L°°(0, T). 


Since Wr is separable, we can boil down a diagonal argument to justify that Jg is relatively compact 
in C°{0,T;W/i — weak — *): we can assume that the extracted subsequence is such that / • 0 dx 

converges uniformly on [0,T] for any 0 S Wr. □ 
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4 Asymptotic analysis of the Vlasov—Poisson—Fokker—Planck 
system 


In this Section we state and prove a Theorem analogous to Theorem 1.2 when the basic equation 
is (VFP), which includes a Fokker-Planck operator, coupled wit h ([P| ). 

For the well-posedness issues of the system (VFP) coupled to (|P|), we refer the reader to nun]. 
The role of the external potential is precisely investigated in m- The associated moment system 
reads 

dtPe + Va, • Je = 0, 
dtJe + • Pe + Pe'^x^e = 

where we still use the notation = J vf^ di;, Pg = J v (S> vfg dz;. As £ —>■ 0, we expect as before 
that Pe —> rie = laA$ext and that the behavior of the current is driven by the Lake Equation with 
friction 

dtv + v- = -y 

Va; • (UeV) = 0. 


(LE/) 


If $ext is quadratic as in (121 (possibly isotropic), the domain 11 is an ellipsoid (possibly a ball) as 
in Section 2.1 and (LE/) becomes the Incompressible Euler system with friction 


dtV + V, • (V 0 V) + = -y 

v. • E = 0. 


(41) 


Eor a more general confining potential 4)ext, we make assumptions hi), h2), HI) and H2) as in 
Section 


2.2 


Since we work with finite charge data, the limit equation (LEj ) holds in H, completed 


with the no flux boundary condition ([^, namely 

V{t,x) ■ v{x] 


on 


= 0 . 


Like in the previous section we associate with V, smooth solution of (LE/), a smooth compactly 

supported extension E defined on [0, T] x such that E • ^{x) = 0. 

dO, 

We shall investigate this asymptotics in the specific case where the “temperature” d — goes 
to 0 as e —>■ 0. In this context, we can derive an analog of Proposition 3.3 that accounts for the 
dissipation mechanisms induced by the Fokker-Planck operator. 

Proposition 4.1 The solution (/e,$e = ^^e+^'ke) of (|VFP|)-(|P|) satisfies the following entropy 
dissipation inequality 


d 

dt 


jfe<ivdx+^ JJ fe 




where we denote 


= JJ \v^/TE + ‘2&e'^v\/%\^dvdx>0. 

Furthermore, the total charge is conserved 


jjf^{t,x,v)dvdx^ jj f^{0,x,v)dvdx 


— m. 


Uniform estimates are not directly included in this statement since the function z z ln(z) changes 
sign. Nevertheless, we can establish such uniform estimates. 
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Corollary 4.2 Wt assume that there exists some (large) A > 1 such that 


y exp(-A$ext)dx < 03. 


(42) 


We suppose also that 0 < e < 1/(8A) and 0 < 9^ < 1. Let x —)■ [0,oo) be a sequence 

of integrable functions that satisfy the following requirements 


du da; = m, 

sup [[ dvdx + d, [[ /f‘|ln(/f‘)| duda; 

8A), 0<6»e<l J J J J 


0<£<1/(8A) 


(43) 






init 1 2 


da;+ J // $e /‘“‘duda; )■ < oo, 


with 


1 




(ne- J fr^dv). 


Let 0 < T < oo and let {fe,^e = be the associated solution of (VFPl -©• Then, 

uniformly for 0 < e < 1/(8A) and 0 < 9g < 1: 

*/ /£{1 + kP + ^£|lii(/£)l) + £~^^efe is bounded in L°°(0,T;L^(M^ x K-^)) , 
a) Vx'^e is bounded m L°°(0,T;L^(M^)), 

Hi) S'e is bounded in L^(0,T). 

Remark 4.3 In any dimension N > \, (|42|) is always true for quadratic potentials. If N = 1 or 


N = 2, hypothesis (42) is satisfied if h2) is, since $ext(a;) + mr(a:) —>■ +oo when |a;| —>■ +oo and 
thatT(x) = —|a:|/2 or — In |a;|/(27r). Therefore, hypothesis needs to be verified only for N > 3. 


Proof. We first observe that hypothesis (42) implies 


J exp(—A$e) da; < oo. 


Indeed, we have $e = F * rie — C"* + d>ext > ‘I’ext — C"*- We write, for h>0, 

fsHf,) < MHf,)\ = /,ln(/,) - 2/,ln(A)(l,- '*<*<1 + lo</e<e-'‘) 


< f,Hf,) + 2hfe + -e-'^/^, 


(44) 


and denote 


Esife) = ^ y/'l'^l^ /eduda;+ ^ J J fe dv dx + ^ J |Va;4'Epda;. 


(45) 


We now use (44) with h{x,v) = \v\'^/{W^) + $e(a;)/(4e0£) to infer 

Oe JJ feMfe)dvdx < 9e JJ /e |ln(/e) | di; da; < 6>e JJ f^ ln{ fe) dv dx + ^E^{fe) 

+ yy exp(-|a)p/(166»E) - $e(a;)/(8e6»E))dudx. 

The last term is equal to 

0£-y exp(—|up/(160e) du y exp(— $e(a;)/(8e0e)) da; < y exp(—|a!p/16) da; J exp(— A$e(a;)) da;. 


31 






thus tends to zero as 0^ —>■ 0 (uniformly for 0 < £ < 1/(8A)). Using the dissipation of the entropy 
given in Proposition |4.1[ we then infer 

+9,JJ fr^\Hfr^)\dvdx > +e,jj /r‘ Hfr^) dvdx 

> E,{f,) + 9,JJ f,\n{f,)dvdx 

> E,if,) + 9,JJ fe\Hfe)\dv dx - |ife(/e) + Ofl,^o(l) 


and the conclusion follows since 0^ < 1. n 


Untill the end of the Section, we shall make hypothesis (42). Since we are dealing with the 
regime 

0 < £ < 1, 0 < 0E < 1, 

the estimates in Proposition 4.2 do not provide L^-weak compactness on the particle distribution 
function and its moments; we still need to work with convergences in spaces of finite measures. The 
first step in the investigation of the asymptotic behavior is summarized in the following claim. 


Lemma 4.4 We make assumptions (42) and (43). Up to a subsequence, we can assume that 
/e converges to f weakly-k in x x M^). Then, converges to = J fdv in 

and in T; — weak-*). Moreover, we can assume that ^ J = 

f vfdv in ./^^([0,T] x M.^), the limit J is divergence-free and supported in [0,T] x U. 


Proof. We follow the arguments of the previous Section. We identify the limit of by coming 
back to the Poisson equation y/eVx ■ Vx'kg = — Pe- The time compactness then appears as a 

consequence of the charge conservation, together with the estimates on the current. We obtain the 
L°°(0,T;L^(K^)) estimate on as in (35). Letting £ go to 0 in the charge conservation equation, 
we obtain dtUe + V^, • J = 0 = Vx ■ J■ Still reproducing the arguments of the previous section, 
based on the conservation of the total charge, we arrive at the following conclusion: 


supp(/) C [0,T] X U X , supp(J) C [0,T] x U. 

Furthermore, J belongs to the set it admits a normal trace J • which actually 

vanishes. □ 

It remains to identify the limit J. As in the case of the pure Vlasov-Poisson equation, the idea 
consists in introducing a suitable functional intended to compare to the expected limit. Let 
-* (0, oo) be a given function such that 


J JT;, dx = m = 


Ug dx 


/(O, X, v) du dx 


and let us set 


Mrfi^{t,x,v) 


1 expf- i^-^(M)n 

( 2 ^ 0,)^/2 ® P 20 , )■ 


A natural candidate to replace the functional would be the relative entropy of /, with respect 

to ng{x)My (t, X, v) associated with the non-negative convex function z ^ z ln(z) — z -|-1, namely 


A In 


fe 


- fe + UgMy^g^ du dx. 


but the first term is clearly meaningless since n, has compact support. Therefore, we introduce 


.XW = ^exp(-—) 


where = 


J exp 


^ejy) 

£ 0 , 


dy, 


(46) 
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and the following modulated functional 

jr/P = dell (^fe In dvdx + ^J dx. (47) 

In fact, is up to the term ^ dx, nothing but the relative entropy of with respect to 

associated with the non-negative convex function G : (0, -l-oo) 9 z !->■ zln(z) — z-|- 1. This 
implies in particular that the integrand in the first integral of ( [47| ) is simply G{fe) — G{NgMy^g ^) — 
G"(/e)(/e — N^My^g^), thus pointwise nonnegative, and vanishes only when = N^My^g^. By 
definition of N,. and My^^ and using the fact Jf f^dvdx = m = ff ^My^g^ dudx in view of our 
normalizations, we infer 


^r!e = JJ fs Mfs) dw dx -f i J j \v - fe dv dx + ^ J 

+ \ I -f ^7Vm6»£ln(27r6»E) - 6»emln 


$e/e du dx 


— 


/e ln(/e) dudx -I- -NmOs ln(27r6'e) - ^^mln 



(48) 


This second expression of justifies the choice we have made for Actually, for our purpose, 
the exact normalization f f ^My^g^ dx dx = m is not necessary, though natural in a modulated 
entropy argument, only the fact that In(yl^My^g^) « —$e(x)/(e0e) — |x — il^p/(2de) is used. This 
is related to the fact that the temperature 0^ is small in the regime we are considering. 

Let us now compare My and My^^ more precisely. As a first step, note that, on the one hand, 

9s ln(27r0e) —^ 0 


when 9s —>■ 0; and on the other hand, that 

|fl| = J exp dy<Zs = J exp dy < J exp (^- A$e(y)) dy < -boo 

if < 1 and e < 1/(8A), thus, as 9s —>■ 0, 


^etnln ( ^ 


0 . 


The inequality (45) implies 

My^s = My^s - JJ fe ln(/e) dx dx - ]^Nm9s \TL{2n9s) + detain 

<^,e-^e JJ /eln(/e)l/,<ldxdx-bOe^o(l) 

< 2My^s + Oe->'0(1)- 

Then, let us compute the time derivative of the modulated entropy My\. We get 

JJ feMfe)dvdx+^ JJ Ix-i^P /edxdx-b iy $e/e dx dx 

+ l J |V,^,|Mx 


(49) 


dt 


= dx dx -b - 


d 

dt 


yjlxp fsdvdxd- i y 4>e/E dx dx-b ^ J |V,j4'£|^dx| 


X- Yfsdvdxd- - 


Vedxdx > . 
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Bearing in mind the computation for proving Proposition |4.1[ we obtain 

d 
dt 


= + J j^.Ydx + ^J Peirpd^rj. 


Reasoning as in the previous section, and by using the moment equations, we are led to 
d 


dt 




+ J {Pef - Je) ■ [dt^ + CV ■ V,)r) dx 
- J oy : {f‘r,e - Va;d;'£ ® V^d'e) dx + i y p^y ■ Va;$e dx, 

using once again that Vx ■ [ngy) = 0 and the no-flux condition Let us set 
^r,e= //'K'y - l^)v^+26»eV„v^|^dxdx > 0. 


We rewrite 


^e = ^r,e- J Pe\y\^dx +2 J y-J.dx. 

Accordingly, we can reorganize terms as follows 

+j{p^y-j^).{dty + {y-Vx)y + y)dx 
- j Dy : (Pr.e - ^ J Pe^ ' dx. 

We can summarize the previous manipulations within the following inequality 

^ + ^r,e <\ j Pe^ ■ V,$e dx + + j Dy : (Pr.e - ® dx, (50) 


where, for any 0 < t < T, 

f- 

Jo 


■ dr = 


{p,y - Je){dry - y ■ v,r + y) dxdr 


tends to 0 as e —>■ 0. We wish to strengthen this result as follows. 


Lemma 4.5 We make assumptions (42) and (43) and suppose that 0^ —>■ 0 as e —>■ 0. We have 

where, for any 0 < t <T, lime_>.o /J dr = 0. 

Proof. We can also reproduce the arguments in the previous section used to estimate 

- / p^y -Vx^edx < - [ Ps^edx. 


For the last term in (50), we have 


J Dr : (Pr,e-V,^e«)V,^e)clx< IlDriloo (^JJ |x-r|2/edxdx + J jV^^el^dx 
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so that, using (49) 


^ JJ fe<^edvdx + i // “ y\'^fedvdx J |V2,4'Epdx + 1 

^ + 0 £-J-o(l)- 


It allows us to conclude by coming back to (50). □ 


Let us now state our main result concerning the Vlasov-Poisson-Fokker-Planck system. We 
recall that we may work either with a quadratic potential $ext (and then the domain 14) is an 
ellipsoid), or with a general potential where hi), h2), HI) and H2) are satisfied. 


Theorem 4.6 If N > 3, we make assumption (42), that is we assume that there exists some 
(large) A > 1 such that 


J exp(-A<i)ext) dec < oo. 


Denote by V the solution, on [0,T ], to the Lake Equation with friction (LEj) with the no-flux 
condition § given by Theorem \A.l\ and consider a smooth extension Y to V. Let /“'* : xK^ —>■ 

[0, oo) be a sequence of integrable functions satisfying 


d?;dx = 


m 


and 


^r,£(/r) ^ 0, 


where is defined in (47). Consider then the associated solutions fe of the Vlasov-Poisson- 

Fokker-Planck equation ( Tv^ -d^. Then, we have, as e ^ 0 and 9^ -y 0, 

i) Pe converges to in C^{0,T; — weak—*); 

a) —>■ 0 uniformly on [0,r]; 


Hi) Je converges to J in x K^), where J\ 


[o,T]xn 


= V, Vx -1 = 0 and J ■ v{x)\g^ = 0. 


Remark 4.7 We have seen (see Remark 4-3) that the integrability assumption J exp(—A$ext) da; < 
00 is automatically satisfied if N = 1, 2 by h2) or for quadratic potentials. When N >3, it is also 
true if 4>ext is convex and tends to +oo at infinity. 

Remark 4.8 One may construct an admissible family of initial conditions following the lines of 
Remark 1.3. In particular, taking G a normalized Gaussian, it is enough to choose 9^ and such 
that 6*£ j /e In —>■ 0, which imposes 0^ In CTe —>■ 0. 


Proof. It is clear that if —)■ 0, then (431 is satisfied. Item i) has already been discussed. 

Applying the Gronwall lemma, we deduce readily that ii) holds from Lemma 4.5 Coming back to 
(49), we infer that JJ \v — 'Pl’^fe drida: tends to 0. Then, we appeal to Lemma 3.5 to conclude that 
J belongs to L°°(0, T; L^(M^)) and that J = UeV. We can also justify some time-compactness as 
in the pure Vlasov-Poisson case. 


A Smooth solutions of the Lake Equations 

Theorem A.l Let Lt be a smooth (dLl of class is enough) bounded open set in j be a 

real constant, s S N such that s > 1 N/2 and rig : t K m 77®^^ such that inf^ Ug > 0. Let 

yinit . ^ ^ divergence free vector field in satisfying the no flux condition I/'"'* • j/ = 0 

on dLl. There exists T > 0 and a unique solution V G L°^{0,T;H^{B{0,R))) of 

r dtV + v-yxV + yxP = -iV, 

\ V, • (ugP) = 0, 
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with the no flux condition §• Moreover, we have 

sup (iiFwiu. + \\dtV{t)\\Hs-i + < c{T) 

0 <t<T '' ' 

for some positive constant C{T) depending on 7 , T, and the initial datum. 


We shall denote V UeV, which is 


Proof. The scheme of proof is exactly the same as in [3 
divergence free. Applying V • (rie .) to the equation, we see that the pressure p satisfies, for any t, 
the elliptic equation 


-Voo ■ xP) = Vx • {UeV ■ VxV) = V 

1 


= y. V, 


= 1^-V., y • V, 


V-Vxl—V 

He 


l<j,k<N 

- Y, djVkdk 

l<j,k<N 



(51) 


where we have used that V is divergence free as well as the identity • (V-VaiZ//)— V• V 2 ,(Va;= 
J2i<j k<N djVkdkUj. We may further impose a suitable Neumann boundary condition for p on 
We recall that for a > N/2, is an algebra. Notice that if y G with s > 1 + iV/2, then the 
right-hand side of (51) is in and 


y-V, 


y-v,(- 

rifi 


E 

l<j,k<N 


d.Vkdk ( ^ 

\ Ufi 


<C\\V\\]j, 




where C depends on infn n-e (which is assumed positive) and the norm of n^. 

Since Ue is in and bounded away from zero and since the boundary is assumed of class 


C®^^, it follows from classical elliptic estimates that (511 endowed with the Neumann condition on 


dV. has a unique solution p G i(l®+^(n) , enjoying the estimate 


lbll^/.+i<c||y||l„ 


(52) 


where C depends on info rig and the i/®"*'^ norm of n^. Assume now that y is a smooth solution 
of (A.l) and let us perform an iJ® estimate. For any a G (N U {0}) with |a| < s, we have 

^ [ |a“yp dx = -2 [ 5“y • ((y • Vx)v) dx-2 [ a“y • o^^Vxp dx -27 [ |a“yp dx 

™ Jq Ja Ja Jq 


Using classical commutator estimates, the Sobolev imbedding i7® C iy^’°° and the iJ®+^ estimate 
we then deduce 

[ | 9 “y |2 dx < - 2 / 9 “y • ((y • Vx)d‘^v) dx + c{\\v\\h^ + \\v\\l. + ||y|| 3 ,.). 

Jci JQ 


(521 on p, we then deduce 
d 


dt , 


We use integration by parts for the first integral (recall that y • ^ = 0 on the boundary), which 
then becomes f^i^x ■ V)\d°‘V\^ dx < Cllylllja. This yields 


^^ |9“yp dx < c(||y|U= + ||y|||,, + ||y||i,.), 

and it follows that, for some Tq > 0 depending only on 7 , Ue and y'"'*, we have ||y||L<»(o,ro;ff'’ )< 
2|!y“"‘||ija. The conclusion of the theorem follows from a suitable viscous approximation where a 
careful treatment of boundary terms is needed, see [33] ■ □ 
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B Construction of an extended divergence—free velocity 

Lemma B.l Let V G L°°{0,T-, {B{0, R),M.^)) be a divergence free vector field in , with 

s > 1 + N/2, satisfying the no flux condition V ■ u = 0 on dB{0,R). There exists a solenoidal 
extension T' of the vector field V defined on the whole space and compactly supported. Namely, 
'f' G L°°(0,T; and it satisfies: 

0 . 

V^-r = 0 in R^. 

Proof. Let us assume TV = 2 or TV = 3. Since -V = 0 in the ball B{0, R), which is convex, there 
exists h G L°°(0,T; ff^~^^(B(0, R),R^)) such that V = V x h. Then, by standard extension results 
(see, e. g., [13 Chapter I: Theorem 2.1 p. 17 & Theorem 8.1 p. 42]), there exists an extension 
h G L°°(0, T; R^)) to h. Considering a cut-off function x G Cf°{R^) such that x(x) = 1 

for X G B{0,3R/2) and denoting 'f = N x (xh), we see that 'f enjoys the desired properties. For 
TV > 4, the construction is similar but involves differential forms. The arguments generalize to the 
case where Q is an ellipsoid. □ 


Remark B.2 In the case of a general potential, Ue is not uniform, and it is possible to construct 
an extension 'f' G {0,T] {R^,R^)) which satisfies: 

0 , 

a) NX ■ (ne'P) = 0 in R^. 

However, it requires further topological hypotheses on LI. Assuming HI, and assuming also that 
K. is connected, and dLl has a finite number of connected components, we may apply \21\ Corollary 
3.2]: UeV is divergence free in the smooth domain fl, hence we can construct a divergence free 
extension ^ : [0,T] x R^ -G R^ to ngV. Using a cut-off function, we may take ^ compactly 
supported in an arbitrary neighborhood of 1C, the latter can be chosen so that A$ext remains > 0. 
Finally, we set P = ^/(A^ext), which is well-defined even when A$ext vanishes. 
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